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VERTICAL TAIL TO THE FREE AERODYNAMIC DAMPING 

OF A MODEL OSCILLATING IN YAW .. 

By John D. Bird, Lewis R. Fisher; 
and Sadie M. Hubbard 


SUMMARY 

The directional damping and stability of a fuselage - vertical-tail 
model oscillating freely in yaw were measured at a Mach number of 0.14 
and compared with the damping and stability obtained by consideration 
of the effects of unsteady lift. The contribution of the vertical tail 
to the damping in yaw of the model agreed within the limits of experi- 
mental accuracy with predictions of the approximate finite-aspect-ratio 
unsteady-lift theories for the range of frequencies tested, and the two- 
dimensional unsteady-lift theory predicted a much greater loss in 
damping with reduction in frequency than was shown by experiment or by 
the approximate finite-aspect-ratio unsteady-lift theories. For fre- 
quencies comparable to those of lateral airplane motions, both the 
unsteady-lift theories and the experimental results indicated that the 
contribution of the vertical tail to the directional stability of the 
model was relatively independent of frequency. 


INTRODUCTION 


The advent of high-speed airplanes of high relative density has 
focused attention on certain problems associated with the dynamic sta- 
bility of airplanes which, because of previous unimportance, have here- 
tofore been neglected. Among these problems is the effect of the perio 
dicity of the airplane motion on the effective values of the various 
stability derivatives. Reference 1 and, more- recently, references 2 
and 3 have indicated the possibility of sizeable effects from this 
problem. 

An appreciable amount of theoretical work on these unsteady-lift 
effects exists at present; however, only a small amount of experimental 
substantiation of the results is. available for the low-frequency range 
of oscillation. As a result a program has been undertaken in the 
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Langley stability tunnel to determine the effects of such variables as 
frequency and amplitude of motion on the contribution of the various 
airplane components to the stability derivatives of present-day air- 
plane configurations. The work reported herein covers that phase of 
the investigation which considers frequency effects on the directional 
damping and stability of a model undergoing a freely damped oscillatory 
yawing motion. The effects of vertical-tail aspect ratio and compres- 
sibility as predicted by the theoretical treatments are discussed in 
relation to the experimental stability characteristics obtained by the 
free-oscillation and by the curved- flow procedures. 


SYMBOLS 


The data are referred to the system of stability axes and are 
presented in the form of standard NACA coefficients of forces and 
moments about the quarter-chord point of the mean aerodynamic chord of 
the normal wing location of the model tested. (See fig. 1.) The 
coefficients and symbols used herein are defined as follows: 

Cl lift coefficient (L/qSy) 

Cj) drag coefficient (D/qSy) 

Cy lateral-force coefficient (Y/qSw) 

C} rolling-moment coefficient (L'/qSybw) 

C m pitching-moment coefficient (M/qSyT^) 

C n yawing-moment coefficient (N/qSwby) 

L lift, pounds 

D drag, pounds 

Y lateral force, pounds 

L' rolling moment, foot-pounds 

M pitching moment, foot-pounds 

N yawing moment, foot-pounds 

q dynamic pressure, pounds per square foot 
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CD 
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N. 

♦ 


mass density of air, slugs per cubic foot 

free-stream velocity, feet per second 

area, square feet 

span, feet 

chord, feet 

mean tail chord, feet 


mean aerodynamic chord, feet 



spanwise distance from the plane of symmetry 
aspect ratio (b^/s) 
angle of attack, degrees 
angle of sideslip, radians 


time, seconds 
yawing-velocity parameter 
angle of yaw, radians 

yawing velocity, radians per second (d\|r/dt) 

yawing acceleration, radians per second per second 
(d 2 +/dt 2 ) 

angular velocity, radians per second 
frequency, cycles per second 

reduced-frequency parameter referred to semichord of 
vertical tail (o>ct/2V) 

mechanical spring constant, foot-pounds per radian 
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t l/2 


m 

P 


Z 

C(k) 
F,G 
P = F 
Q = H 
a 

Z,M 

A,B 

C nr = 


yawing moment of inertia, foot-pound- second^ 
time to damp to one-half amplitude, seconds 

damping constant 
period of oscillation, seconds 

nondimens ional tail length referred to semichord of 

/ Z 

vertical tail 

\ c t/2 

distance from origin of axes to midchord point of vertical 
tail, feet 


= F + iG 


+ iG 
+ iJ 


^£n 
a rb 
^ 2V 


circulation functions used by Theodorsen (reference 6) 

finite-span functions used by Biot and Boehnlein 
(reference 9) 

finite- span correction used by Reissner (reference 8) 

i 

compressible-flow functions used in place of F and G 
(references 10 and ll) 

functions of k defined in text 
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Subscripts: 
w wing 

t vertical tail 

0 amplitude 

f due to friction 


APPARATUS 

Model and Oscillating Strut 


The model used in this investigation and shown in figure 2 had 
separable wing, fuselage, and tail surfaces and was constructed almost 
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entirely of balsa wood in order to minimize the mass and make the 
natural frequency of the model on its mounting as high as possible. 

The configuration is representative of that of several present-day high- 
speed airplanes. The aspect ratio A and mean chord ct of the 
vertical tail were 2 and 4.9 inches, respectively. 

The model was mounted at the quarter chord of the wing mean aero- 
dynamic chord on a 1- inch-diameter rod which rotated with respect to a 
circular hollow strut and which was connected to it by means of flexure 
pivots made of Swedish steel. Flexure pivots were used rather than 
bearings in order to minimize the friction in the system. The pivots 
were alined with the center of the rod on which the model was mounted. 
The entire assembly, which is shown in figures 3 and 4 , was fastened to 
the ceiling of the tunnel test section and covered by a streamlined 
fairing. 

The inner rod or oscillating strut, which was permitted to rotate 
and upon which the model was mounted, extended through the ceiling of 
the tunnel test section and was capped by a horizontal crosspiece. -The 
period of oscillation of the model was controlled by varying the moment 
of inertia about the vertical axis of the system. This variation was 
accomplished by using crosspieces of different lengths and by clamping 
weights at the ends of the crosspieces. 


Recording Device 

A continuous record of the displacement in yaw of the model was 
made on photographic film by means of an optical recording system. A 
mirror clamped to the section of the oscillating strut which extended 
outside the tunnel reflected a beam of light from a point light source 
within the recorder. The beam was focused on a roll of photographic 
film by means of an 8- inch-focal-length double -concave lens. The film 
roll was driven at constant speed by gearing it to an electric motor. A 
10-cycle -per- second neon timing light incorporated in the recorder 
exposed timing lines on the film in order that time as well as model 
displacement could be measured accurately on the film record. 


TESTS 

Oscillation Tests 


The free-damping tests were conducted with the model in an inverted 
position, in order to minimize strut interference effects, and at zero 
angle of attack with respect to the fuselage reference line. The tests 
were conducted on the combination of fuselage and tail surface and on 
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the fuselage alone "because preliminary tests had shown that the presence 
of the wing had no apparent effect on the damping at zero angle of 
attack. The range of moments of inertia and frequencies covered are 
shown together with the test results in table I. Results of runs made 
with wind on and off are given in order to illustrate the magnitude 
of the tares applied. 

When tests were made of the fuselage alone, a torque rod was 
employed to provide a stable restoring moment to the fuselage which 
was equal to the aerodynamic moment provided by the vertical tail when 
attached. The torque rod was clamped between the extended oscillating 
strut and a building structural member as shown in figure 3. The 
spring constant of the system was 6.8 foot-pounds per radian for the 
flexure pivots alone and 32.8 foot-pounds per radian for the combina- 
tion of the flexure pivots and the torque rod. 


Force Tests 

The steady-state stability characteristics of the model were deter- 
mined through an angle-of-attack range from 0° to 8° for the fuselage 
alone, the fuselage and tail surfaces, and the complete configuration 
by means of wind-tunnel force tests. The model was mounted on a single- 
strut support and all forces and moments were measured by a six-component 
balance system. The steady-state sideslipping derivatives were deter- 
mined through a sideslip range of ±5°, and the yawing derivatives, at 
values of the yawing- velocity parameter rb/2V of 0, -0.0291, -O.0616, 
and -0.08l0 by standard stability- tunnel curved-flow testing procedures. 

All tests were conducted in the 6- by 6-foot test section of the 
Langley stability tunnel at a dynamic pressure of 24.9 pounds per square 
foot, a Mach number of about 0.l4, and a Reynolds number of 442,000 
based on the wing mean aerodynamic chord. Standard jet-boundary correc- 
tions have been applied to the angle-of-attack, the drag-coefficient, 
and the rolling-moment-coefficient data obtained from the force tests. 

In addition, the steady-state yawing data have been corrected for the 
effects of the lateral static -pres sure gradient peculiar to the curved- 
flow testing procedure. Jet-boundary corrections were not applied to 
the oscillation tests because of their insignificance for the test 
conditions involved. See reference 4. 


REDUCTION OF OSCILLATION TEST DATA 


From the continuous film record taken of the motion of the model 
after an initial displacement, the amplitudes of the successive cycles 
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were measured and plotted to a logarithmic scale against time. Inasmuch 
as the damping is of a logarithmic nature, the resulting plot is a 
straight line from which may be read the time for the motion to damp to 
one -half amplitude ^ 1 / 2 ) • 

For a system having a single degree of freedom in yaw such that 
for the stability-axis system the angle of sideslip p is the negative 
of the angle of yaw. \|r, the equation of motion can be expressed as 


. ( x z - »y)V " ( N * “ n p + . N * f )* >• (S - N p)* = 0 . d) 

where and refer to the mechanical friction of the system and 

the flexure-pivot spring constant, respectively. The solution of 
equation (l) can he written as 


where 


\|f = -e _irt ('A sin 2itft + B cos 2nft) 


2jtf =1 


/% - Np 


- + N 

if 

y i z - 

_ 2 I 

pz - N V) 

l 
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m = -- 


- “e + \ 

i 1 * 


( 2 ) 




The damping constant may he expressed in terms of the time to damp to 
one-half amplitude as 


m 


-0.693 

t l/2". 


The directional stability and damping of the system may then he 
expressed in nondimensional form in a manner similar to that of 
reference 5 as 





( 3 ) 
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-2.772(I Z - H$)V 
qS w b w 2 




wind off 


00 


Iii order to find the damping due to friction in the system, the damping 
of the model is determined at zero wind velocity with an equivalent 
mass in place of the vertical tail and a torque rod of sufficient 
strength to produce the proper frequency, inserted in the system. The 
term N^r is interpreted in terms of a factor added to C n ^ because 

this particular group of free-damping tests afforded no convenient 
means for its separation. For the frequencies employed in these tests 

/' b w\ 2 


the factor C 


f k 


V, 


should he small in comparison with C 


np- 


The damping derivatives were obtained from data similar to the 
sample records of model motion and plots of logarithm of amplitude 
against time given as -figure 5- A small residual yawing motion of the 
model attributable to turbulence in the tunnel air stream necessitated 
determining the time to damp to one-half amplitude by the slope of the 
amplitude envelope at the largest angles of yaw in order to minimize 
the error in determining the damping derivatives. This source is 
believed to have resulted in no more than a 5-percent uncertainty in 
the results. 


ANALYTICAL CONSIDERATIONS 


Various unsteady-lift theories were employed to calculate the 
contributions of the vertical tail to the directional stability and 
damping of the model for comparison with the experimental results. 
These theories included the incompressible cases covered by Theodorsen 
for two dimensions and by Jones (reference 7) > Reissner (reference 8), 
and Biot and Boehnlein (reference 9) with various degrees of approxi- 
mation for finite span. The two-dimensional work including the sub- 
sonic effects of compressibility as given by Possio and Frazer and 
Skan is also employed. A short summary indicating procedures for 
obtaining the necessary derivatives from these theories is included 
here for completeness. ’ .... 


Two-Dimensional Derivatives 

The work of Theodorsen on the derivation of the expressions for 
the unsteady lift and moment of a two-dimensional surface undergoing 
sinusoidal oscillations in an incompressible fluid is given in 
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reference 6. From this source, the moment per unit span about the air- 
plane center of gravity of a vertical tail mounted on an airplane which 
is performing sideslipping oscillations can he shown to he 


NO) = -irqct 2 


ac t 3 


W 


(a + |)c(k)3 


(5) 


A similar expression for an airplane performing yawing oscillations is 


N(t) = -% qct 3 * 




( 6 ) 


An appendix to this paper gives a discussion of the conversion of these 
results from reference 6. 

The Theodorsen function C(k) is a complex circulation function of 
the reduced frequency and in terms of its real and imaginary parts is 
written 


C(k) = F + iG 

The functions F and G may he computed from expressions given in 
reference 6 in terms of Bessel functions of the first and second kinds 
of orders zero and one and of argument k. In stability practice equa- 
tions (5) and (6) are usually interpreted in terms of derivatives as 
follows: 


N(P) 


qs w b 2 • 

Cn p qS w h w 3 + c n - — Zf- 3 


NO) 
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qSwhw^ 

2V 


♦ + c T 


qSyhw 3 

kv 2 
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3 = 3 0 e 


icot 


i|r = \|r 0 e 


iojt 


If the sinusoidal motions 
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are employed for a substitution 
real and imaginary terms may be 

C nr , and C n £ determined. The 
follows: 


into equations ( 5 ); ( 6 ) and ( 7 ), the 
equated and the derivatives C n p, c n^ 

results of such an analysis would be as 



The velocity terms contain only the F term of the Theodorsen function 
whereas the acceleration terms contain the G term divided by k. The 
corresponding side-force derivatives are derived in the appendix. For 
the free-oscillation tests conducted for this paper, the motion is such 
that the angle of sideslip is the negative of the angle of yaw. The 
damping derivative measured therefore is the difference between the 
derivatives C n ^ and C n p determined from equation ( 8 ) . This fact 

was illustrated previously with respect to the reduction of the oscil- 
lation test data. In. a like manner, the directional stability parameter 
measured can be shown to be a combination of the C n p and C n ^ terms, 

specifically 

% * 

Combination of the results of equations ( 8 ) to obtain the factors 
calculated for comparison with the oscillation tests gives 
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where 



(9) 


Finite-Span Derivatives 

The aerodynamic -span effect is considered in reference 7 for wings 
of aspect ratios of 6 and 3 by correcting the aerodynamic inertia and 
the angle of attack of the infinite-span surface. An approximation 
employed in this analysis concerning motions of long wave length may 
make the results subject to question for values of k near zero. The 
reference presents expressions for calculating the finite-span values 
of the Theodorsen function C(k). The use of these F and G circu- 
lation functions in equations (8) and (9) results in values of the 
various derivatives considered for aspect ratios of 6 and 3* 

The aerodynamic- span effect is considered for wings of arbitrary 
aspect ratio in reference 8. The three-dimensional effect of the 
finite span may be obtained by adding a correction term a to the 
basic Theodorsen two-dimensional function. The finite-span func- 
tion C(k) + a is employed in a rather lengthy and tedious computation 
in order to determine the components of the unsteady lift for each of 
five points along the span of the oscillating surface. The various 
stability derivatives may then be determined by a graphical integration 
of the in-phase and out-of-phase components. An approximation made in 
this analysis is of such a nature as to make the method primarily 
suited to the higher aspect ratios. A practical lower limit of appli- 
cability is not known. 

A one-point approximation to the span effect for the calculation 
of the unsteady lift may be employed rather than the five-point proce- 
dure. Comparative calculations made for the present investigation 
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have shown the one-point approximation to yield results for low aspect 
ratios which are essentially the same as those calculated "by the 
longer procedure. 

Biot and Boehnlein (reference 9) considered the aerodynamic- span 
effect on unsteady lift by an approach that yields a closed form for a 
one-point approximation to the span effect. The finite-span circula- 
tion functions used are 


P = F + iG 


Q = H + iJ 


For infinite aspect ratio 


P = Q = C(k) 


where C(k) is the two-dimensional Theodorsen function. The func- 
tions P and Q may be calculated for the argument k by methods 
given in reference 9, and equations (8) and (9) are again used to 
calculate the various stability derivatives. The factors A and B _ 
which appear in equations (9) must, however, be written in terms of P 
and Q as follows: 

kJ + [ a ^ + 


B = (2a - l)akF + (a - |)kH - (a - -ijk - 2aG - J 


5) 


1st + 2aF + H 


A = (2a - l)akG + 


Effect of Compressibility 

Subsonic compressibility effects on the unsteady damping in yaw 
and directional stability, of a two-dimensional oscillating lifting 
surface can be considered through the use of reference 10 which is an 
extension of the calculations of Possio (reference 11). The infinite- 
span compressible-flow values for the stability derivatives can be 
obtained by the use of equations (8) and (9) for Mach numbers 0, 0.5, 
and 0.7- In the notation of Frazer and Skan (reference 10), however, 


A = 

.1(2 

2V2 

; + ^ + z 3 

-I z i) 

+ M 

B = 

2V2 

j + + Zij. 

-i z 2) 

+ M, 

i 
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The correspondence between the notation of Theodor sen and that of 
Frazer and Skan is shown in reference 12. The values of the Z and 
M functions are tabulated by Frazer and Skan for Mach numbers 0 , 0 . 5 , 
and 0.7 and a significant range of reduced frequencies. Reissner has 
recently extended his work on oscillating wings of finite span to 
include, to the same degree of approximation as for the incompressible 
case, the subsonic effects of compressibility. (See reference 13 . J 


RESULTS AND DISCUSSION 
Presentation of Results 


The experimental values of the longitudinal, directional, and 
yawing characteristics .of the complete model, the model with wing 
removed and the fuselage alone as determined by force tests are given 
in figure 6 . The experimental values of the damping in yaw of the 
fuselage and tail assembly (designated FT) and the fuselage alone 
(designated F) obtained by the freely damped oscillation technique for 
several frequencies of oscillation are shown in figures 7 and 8 . The 
frequency and aspect-ratio effects on the damping in yaw contributed by 
the vertical tail as calculated from the unsteady-lift theories dis- 
cussed previously are also shown in these figures. Figure 9 gives a 
comparison of the relative importance of Cn r and C n p on the damping 

in yaw of the model. Figures 10 to 12 are similar figures summarizing 
frequency and aspect-ratio effects on directional stability. Figures 13 
and l4 give the computed effects of compressibility on the vertical- 
tail contribution to the damping in yaw and directional stability of 
the test model. 


Discussion 

The variation of the experimentally determined damping in yaw of 
the test model with frequency is in reasonably good agreement for the 
range covered by the tests with the various approximate theoretical 
treatments which consider the effect of finite span. (See figs; 7 and 8 .) 
The contribution of the fuselage to the damping in yaw is small for the 
range of test frequencies. A two-dimensional approximation to the 
damping as is frequently employed for flutter work appears to be 
remarkably good for k > 0.1 and for aspect ratios of 6 and greater 
(fig. 8 -). A large effect of aspect ratio, however, is shown for aspect 
ratios below 6 for the entire frequency range and for all aspect ratios 
for frequencies less than k = 0.1. The two-dimensional result indi- 
cates a loss in damping at low frequencies much greater than the finite- 
span result and, in fact, shows a reversal in sign of the damping at 
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sufficiently low frequencies. A number of papers dealing with unsteady 
lift have noted this characteristic of the two-dimensional aerodynamic 
forces. (See references 1 to 3 and 1^.) 


An examination of the values of the derivatives C n and C n * as 

r 

computed from reference 6 for infinite aspect ratio and' reference 7 for 
aspect ratios of 3 and 6 (fig. 9) indicates that the derivative C n ^ is 

critically dependent on reduced frequency for the low range of reduced 


frequencies and infinite aspect ratio. The derivative C n is only 


slightly dependent bn the reduced frequency for infinite aspect ratio 
and even less so for the lower aspect ratios. The dependence of C n p 

on reduced frequency is responsible for the decrease in damping shown 
for the infinite-aspect-ratio case as the frequency is reduced (figs. 7 
and 8) . 


The aspect-ratio-2 results ' shown in figures 7 and 10 were obtained 
from the theory of Jones (reference 7) for aspect ratio 3 by use of a 
correction factor which was the ratio of the lifting- line- theory lift- 
curve slopes for aspect ratios of 2 and 3* This correction is only 
approximate and was made in order to place the results on a comparable 


basis. The variation of 


with reduced frequency, shown in 


figure 7> indicates no appreciable effect of frequency to zero frequency. 
In fact, the result of Jones for aspect ratio 3 reduces to a finite 
value at k = 0 given by 


jt F t \ £t , 

'aW, Sw K 


1.2a^ - 0.370a + 0.515) 


which is not the case for the two-dimensional result that approaches 
infinity -as k approaches zero. This result is not expected to he 
very accurate close to k = 0, however, because of the nature of the 
assumptions made in its derivation as mentioned previously. ■.* 

When the finite-span effects are considered, it appears that these 
effects. are especially significant with regard to the damping in yaw of 
this model for aspect ratios smaller than 6. (See fig, 8.) The theory 
also indicates a significant decrease in the damping at low frequencies 
when the aspect ratio was increased from 20 to-oo. 

The experimental values of damping in yaw obtained by free oscil- 
lation agree closely with the values of C n obtained by the standard 

curved-flow testing procedure used in the Langley stability tunnel. The 
curved-flow technique results in the k = 0 value of C n . No unsteady 
effects, however, are obtained by this technique. 1 
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The unsteady-lift theories indicate that, at frequencies of oscil- 
lation from k = 0 to k = 0.1, the frequency has little effect on the 
directional stability contributed by the vertical tail but that high 
frequencies of oscillation can cause the directional stability con- 
tributed by the vertical tail to decrease greatly and eventually to 
change sign (figs. 10 and ll) . Changes in aspect ratio, for low aspect 
ratios, have a large effect on the directional stability; for high 
aspect ratios, changes in aspect ratio are less significant. The 
experimental values of directional stability obtained from the free- 
damping tests remained constant for the range of frequencies investi- 
gated, although the magnitudes are somewhat smaller than those predicted 
for the tail contribution by the calculation methods for aspect ratio 2. 
This result is not unexpected because all the finite-span unsteady-lift 
theories employed maintain some resemblance to the Prandtl lifting-line 
theory which, of course, yields values of the lift-curve slope for low 
aspect ratios which are unrealistically high. The values of directional 
stability obtained by both the free -oscillation and curved- flow proce- 
dures were about the same magnitude. 


Figure 12 gives a comparison of the contribution . of a vertical tail 
to the stability derivatives C no and as computed from refer- 


' n 3 


ences 6 and 7- These derivatives were combined in the proportion 


M 2 


' n (3 


+ C n .k |— 1 


: t/ 


for comparison with the results of the oscillation tests. Although C n . 

V 2 


is of appreciable magnitude itself, the factor C 




is very 


small relative to C n ^ for the range of frequencies around k = 0.01 

and thus should not greatly affect the directional stability of the * 
usual airplane. - 


The computed two-dimensional Mach number effects on the damping in 
yaw and directional stability of the tail surface of the test model are 
shown in figures 13 and 1^. These results were included to illustrate 
the powerful . effect of compressibility on the damping contributed by 
the tail surface. A complete coverage of one-degree-of-freedom flutter 
by use of these aerodynamic forces which may be obtained from refer- 
ences 10 to 12 is given in reference l4. Compressibility effects 
decrease the tail-surface damping at low frequencies and increase the 
damping at higher frequencies. The directional stability increases at 
all frequencies with increase in Mach number with the largest increases 
occurring in the high-frequency range of oscillations. Reference 15. 
which presents the two-dimensional oscillating aerodynamic forces at 
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sonic speed indicates larger deviations from the results for incompress- 
ible flow than for the subsonic results given here. 


CONCLUSIONS 


The following conclusions were drawn from the results of an inves- 
tigation made at a Mach number of 0.1 1 * of the effects of frequency and 
aspect ratio on the directional damping and stability of a fuselage and 
vertical tail oscillating freely in yaw: 

1. The contribution of the vertical tail to the damping in yaw of 
the model agreed within the limits of experimental accuracy with predic- 
tions of the approximate finite-aspect- ratio unsteady-lift theories for 
the range of frequencies tested. 

2. Two-dimensional unsteady-lift theory predicted a much greater 
loss in damping with reduction in frequency than was shown by experi- 
ment or by the approximate finite-aspect-ratio unsteady-lift theories. 

3. Both the unsteady-lift theories and. the experimental results 
indicated that the contribution of the vertical tail to the directional 
stability of the model was relatively independent of frequency for 
frequencies comparable to those of lateral airplane motions. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., November 2 6 , 1951 
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APPENDIX 


CONVERSION OF UNSTEADY FORCES AND MOMENTS FROM 
FLUTTER TO STABILITY AXES 


The conversion of the contribution of the vertical tail to the 
yawing moment of an airplane in terms of unsteady-lift functions and 
the Theodorsen flutter axes to the stability axes which are fixed in 
the body is included here as a matter of general interest. Reference 16 
gives considerable attention to such axis conversions. The following 
sketches indicate the paths and attitudes of airplanes performing 
sinusoidal sideslipping and yawing dscillations within the limits of 
the stability-axis system: 
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The Theodorsen coordinates h and a and the stability-axis coordi- 
nates P and t are shown with respect to the paths followed by the 
airplanes in sideslipping and yawing flight. 

In sideslipping flight 



a = 0 


and in yawing flight 





a 


From reference 
of an airplane 

N(a,h) 


6 the contribution of the vertical tail to the yawing moment 
in terms of the flutter coordinates is 



If this moment is for a sideslipping motion or a yawing motion in the 
sense given by stability axes fixed in the body and the following terms 
are substituted for the sideslipping case 


a = 0 
a = 0 
a = 0 

h = -Vp 

• ♦ • 

h = -Vp 
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and for the yawing case 


fi = -Vt 

h = -V\Jf 


a = 

• • 

a = f 

• • •• 

a = t 

the following results are obtained for the sideslipping case 

pc, 2 /aitc,V \ pVc 2 tt . . 

N(P) ^(—2^ P ) + — |~( a + i)c(k)(-VP) 

and for the yawing case 

.(♦) - # * ♦ ♦ * 2 )* + ^ ♦] * 


pV c . 


n(a + - ajifrj 


which reduce for the sideslipping case to 


NO) = -*qc t 2 ( + (a + |)c(k)pj 


and for the yawing case to 


»<♦> =-f ^t 3 -! - cw(i - * 2 )]l + t(§ + *% 


Similar expressions may be obtained for the side force in pure side- 
slipping and yawing flight. Reference 6 gives the side force of the 
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vertical tail in terms of the flutter axes as 

2 , - x 
P c t / .. flc+aaX 1- 

Y(a,h) - — (vna + ith - — ^ — I - itpVc^.C(k) Va + h 

which, when the following terms are substituted for the sideslipping 
case 

h = -VP 
h = -VP 



and for the yawing case 


a = 0 
a = 0 
a = 0 

h = -V* 
h = -V* 


a = i|r 

• m 

a = )| f 

• • • • 

a = i|r 

becomes for pure sideslipping and yawing flight 

Y(P) = qc t f[c t (|)+ ^(k)p| 

and 


,Y(t) 
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Substituting the motions, in the manner done in the text, for the 
yawing moment, 


into these expressions and 
derivatives which are 


Y(P) 


yU) 


3 = 3 0 e 


ioyt 

inrt 


= \|f e 
0 

into the expressions in terms of stability 
• ^w^w 

= + Cy.3 - 2 ^— 


= c Y r ^ 


qS w bw 


2V 



4v 2 


gives the following relationships for the tail contribution to the side- 
force derivatives when the real and imaginary parts of these equations 
are equated: 

Page 22: Last four equations should be: 
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Figure 1.- System of stability axes. Arrows indicate positive forces, 
moments, and angular displacements. Yaw reference generally chosen 
to coincide with initial relative wind. 
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Figure 3*- Diagram of test equipment. 
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Figure k.- Test model mounted on free-oscillation strut in Langley 

stability tunnel. 



Increasing time 
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Figure 10.- The effect of frequency on the directional stability for i 
low range of the reduced frequency. Comparison of several methods 
for calculating frequency effect on C n . 
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method of Biot for a range of reduced frequencies. 



(b) Range of reduced frequencies. 

Figure 13.- The effect of compressibility on the damping in yaw 






